Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species.

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Mitochondrial genomes of the freshwater monogonont rotifer Brachionus fernandoi and of two additional B. calyciflorus sensu stricto lineages from Germany and the USA (Rotifera, Brachionidae).

The Brachionus calyciflorus species complex was recently subdivided into four species, but genetic resources to resolve phylogenetic relationships within this complex are still lacking. We provide two complete mitochondrial (mt) genomes from B. calyciflorus sensu stricto (Germany, USA) and the mt coding sequences (cds) from a German B. fernandoi. Phylogenetic analysis placed our B. calyciflorus sensu stricto strains close to the published genomes of B. calyciflorus, forming the putative sister species to B. fernandoi. Global representatives of B. calyciflorus sensu stricto (i.e. Europe, USA, and China) are genetically closer related to each other than to B. fernandoi (average pairwise nucleotide diversity 0.079 intraspecific vs. 0.254 interspecific).